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Abstract 

We review and update the effects of physics beyond the standard model on CP asymme- 
tries in B decays. These asymmetries can be significantly altered if there are important 
new-physics contributions to B®-B® mixing. This same new physics will therefore also 
contribute to rare, flavor-changing B decays. Through a study of such decays, we show 
that it is possible to partially distinguish the different models of new physics. 
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1. Introduction 



Within the standard model (SM), CP violation is due to nonzero complex phases 
in the Cabibbo-Kobayashi-Maskawa (CKM) quark mixing matrix. In the Wolfenstein 
parametrization [1], only the elements V u b and Vtd have non- negligible phases: 



where A = 0.22 is the Cabibbo angle. The phase information of the CKM matrix can be 
displayed elegantly using the so-called unitarity triangle (Fig. 1), which follows from the 
orthogonality of the first and third columns. CP violation is indicated by a nonzero area 
of the unitarity triangle; to date, the only evidence for CP violation comes from \e K \ in 
the kaon system. 



Fig. 1: The unitarity triangle. The angles a, (3 and 7 can be measured via CP violation in the B 
system. 

At present, constraints on the unitarity triangle come from a variety of sources [2]. 
The sides of the triangle can be probed directly - (K^/V^l in charmless B decays, and 
I^WK&I through B^-B^ mixing. The three angles, a, (3 and 7, are constrained by the above 
measurements, as well as those of \e K \ and B decays to charmed mesons (|V c b|). However, 
with the exception of \V c b\, in all cases there are large theoretical hadronic uncertainties 
in the extraction of the CKM matrix parameters from such measurements. As such, our 
current knowledge of the unitarity triangle is rather poor. 
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In the coming years, we will be able to precisely determine the unitarity triangle, 
and hence test the SM explanation of CP violation. The key measurements involve CP- 
violating asymmetries in B decays. Through such measurements, the weak phases a, (3 
and 7 can be extracted with no hadronic uncertainty [3] , and then compared with the SM 
predictions. 

It is useful at this point to review how the CP angles are probed in CP asymmetries. 
Most asymmetries of interest measure mixing-induced indirect CP violation, which comes 
about through the interference of the two amplitudes B — > / and B — > B — > / '. In order 
to cleanly extract the weak phases from these CP asymmetries, two conditions must be 
met. First, in the neutral B system, one must have r 12 <C M 12 . This relation holds within 
the SM, where r 12 /M 12 ~ 37r m 2 /m 2 < 10~ 2 . Second, the direct decay B — > / must be 
dominated by a single weak amplitude. If this is not the case, then one may have direct 
CP violation, which involves unknown strong phases. In fact, in the SM most B decays 
which are useful for CP asymmetries have more than one weak amplitude - in addition to 
the tree-level contribution, one may also have penguin diagrams [4]. However, for the cases 
of interest, the penguin contamination is either unimportant or can be eliminated using 
isospin [5] and other considerations. We refer to Ref. [3] for a more complete discussion of 
these issues. 

Assuming that the above two conditions are met, the CP asymmetry measures ImA, 
where A is a pure phase: 



The three pieces are defined as follows [6] . X is the weak phase of the direct B — > / decay 
amplitude. For example, X^ Uu ^ = V u bV* d in the SM. Y is the phase of B°-B° mixing: 
e.g. for B®-B® mixing in the SM, Yd = V t * b V t d- Finally, Z is the phase of K-K mixing, 
which is important only if the final state / contains a neutral kaon. In the SM, assuming 
that K-K mixing is dominated by box diagrams with virtual c quarks, Z = V c dV* s , which 
is real to a good approximation in the Wolfenstein parametrization. 

From the above equation, it is straightforward to establish which CP angles are mea- 
sured in different CP asymmetries. For example, the CP asymmetries in B d — > n + n~ and 
Bd — > ^K s probe sin 2a and sin 2/3, respectively. And the angle 7 can be extracted from 
the CP asymmetry in B s — > DfK^ [7] (the function in this case is sin 2 7). Another way 
of measuring 7, which doesn't involve mixing-induced CP violation, is via the asymmetry 
in B^ — > D CP K ± [8]. In all cases, the CP phases can be obtained with no hadronic 
uncertainty. Of course, there are many other CP asymmetries which can be used to obtain 
the angles a, (3 and 7. (For certain decays {e.g. B s — > ^(fi) CP asymmetries probe very 
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small angles of other unitarity triangles, which are almost flat [9].) 

Once these angles are measured, it will be possible to test the SM by comparing the 
measured values with the SM predictions, as well as with the angles expected from inde- 
pendent measurements of the sides of the unitarity triangle. (As we will see below, these 
two comparisons are not necessarily equivalent.) As mentioned above, the SM predictions 
are not very precise at present. Even so, the experimental data do somewhat constrain the 
CP angles [2]. For example, sin 2/3 must be between 0.32 and 0.94 at 95% c.l. In addition, 
the predictions for a, (3 and 7 are correlated, since there is only a single complex phase 
in the CKM matrix and the three angles must add up to 180°. A special correlation was 
shown to exist between small values of sin 2/3 and large values of sin 2a [10], and an almost 
linear correlation was found between a and 7 [11]. 

All this presupposes that the SM is the complete description of the weak interactions 
and CP violation. However, it is widely accepted that there must be physics beyond the 
SM. There are a number of ways in which new physics can manifest itself through the 
measurement of CP asymmetries: 

• (1): The relation a + /3 + 7 = 7ris violated. 

• (2): Although a + (3 + 7 = 7r, one finds values for the CP phases which are outside of 
the SM predictions. 

• (3): The CP angles measured are consistent with the SM predictions, and add up to 
180°, but are inconsistent with the measurements of the sides of the unitarity triangle. 

In any of these cases, it is only natural to then ask what type of new physics could be 
responsible. (A special possibility is the so-called superweak-type model [12], in which the 
CKM phase vanishes (i.e. 7 = 0), and the CP asymmetries in Bd — > n + n~ and Ba — > ^>K S 
are equal in magnitude. It is possible, but not necessary, that these asymmetries vanish 
in such models, in which case the unitarity triangle becomes a straight line.) 

A first step in answering this question was taken in Ref. [13], which examined the 
effects of a variety of models of new physics on the CP asymmetries in B decays. The 
authors of Ref. [13] concentrated on items (1) and (2) above. Their conclusion was that 
the predictions of the SM can be considerably altered in many of these models. 

The reasoning goes as follows. First, it is very difficult to significantly change the 
relation T12 <C -^12, even in the presence of new physics. Second, in most models of new 
physics, there are no new tree-level contributions to B decays. Therefore CP asymmetries 
continue to measure a well-defined CP phase, as in Eq. (2). Of the three pieces in Eq. (2), 
only Y is likely to be significantly altered by new physics. This is because (i) although 
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Z may be modified in the presence of new physics, only in extremely contrived models 
can arg(Z) be changed, and (ii) X can be affected only if there are new amplitudes which 
can compete with the VF-mediated tree-level decay, and there are very few models of new 
physics in which this occurs [14]. Thus, the principal way that the SM predictions for CP 
asymmetries can be significantly modified is if there are sizeable new-physics contributions 
to B°-B° mixing with phases different than in the SM. It is therefore straightforward to 
establish, model by model, which types of new physics can do this. Examples of models of 
new physics which can significantly affect the CP asymmetries include Z-mediated flavor- 
changing neutral currents, four generations, nonminimal supersymmetric models, etc. In 
some of these models, \€ K | is also likely to obtain sizeable new-physics contributions. 

One point which was not emphasized in Ref. [13] is the third way of detecting new 
physics (item (3) above). That is, even if the new-physics contributions to B mixing have 
the same phase as in the SM, their presence can still be detected. This is because, although 
the CP phases a, (3 and 7 are unchanged from their SM values, the new physics affects 
one of the sides of the unitarity triangle, namely the extraction of (V^/V^l from B d -B^ 
mixing. Thus, the measurements of the angles and those of the sides will be inconsistent 
with one another, indicating the presence of new physics. There are several models in 
which this may occur - two-Higgs-doublet models with natural flavor conservation, minimal 
supersymmetric models, etc. 

The question of new physics and CP asymmetries in B decays has also been discussed 
in Ref. [15]. This paper focussed on how new physics can affect various relations among 
CP asymmetries in the SM. One of the points made, which is of particular interest for our 
purposes, is the following. Suppose that a (which stands for n — f3 — 7) and f3 are measured 
via CP asymmetries involving B% decays, and that 7 is obtained through B° s decays. In 
this case, if the phase in B° s -B® mixing is identical to that of the SM, then the relation 
a + (3 + 7 = 7T will hold, regardless of whether there is new physics in B° d -B° d or B° s -B° s 
mixing. In other words, any new-physics effects in B d -B d mixing cancel in the sum of a 
and p. 

This has important experimental implications. The angles a and f3 will probably 
be measured through CP asymmetries in B d — > n + n~ and B d — > ^>K S , respectively. If 
the angle 7 is measured through the CP asymmetry in B s — > DfK^, then one might 
find that the three CP angles do not add up to 180°, if there is new physics in B®-B® 
mixing. However, if the angle 7 is obtained via B^ 1 — > D CP K ± , then, unless there are new 
contributions to B decays, one must find ct + /3 + 7 = n. This underlines the importance of 
measuring 7 (as well as a and (3) in a variety of independent ways. This also demonstrates 
that it will be crucial to search for new-physics effects in all three ways - if only one of the 
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methods (l)-(3) is used, one might miss the presence of physics beyond the SM. 

This rather lengthy introduction summarizes previous work on new physics and CP 
asymmetries in the B system. However, these analyses only partially address the issue. 
Suppose that the CP angles a, f3 and 7 are measured, and it is found that, in fact, the 
presence of new physics is indicated. Ref. [15] presents some tests to determine where the 
new physics might be found {e.g. B^-B® mixing, B° s -B° mixing, etc.) Ref. [13] identifies 
which models of new physics could be involved. However, neither of these references tells 
us how to distinguish among these various models. It is this question which we address in 
this paper. 

As argued above, the new physics can affect the CP asymmetries mainly through 
its contributions to B^-B^ or B®-B® mixing, which are flavor-changing processes. This 
same new physics will therefore also affect rare flavor-changing decays, such as b — > sX 
or b — » dX. (In this paper we generically refer to such processes as "penguin" decays.) 
This is the key point. As we will show, some models of new physics can be distinguished 
by their contributions to these rare processes. In fact, for certain models, the new-physics 
parameter space leading to large contributions to B-B mixing also predicts large deviations 
from the SM predictions for certain penguin decays. Conversely, if no deviation from the 
SM is found, this would so constrain the parameters of the new physics as to render its 
effects in B-B mixing, and hence the CP asymmetries, unimportant. It is an experimental 
question whether or not measurements of the rates for such penguin decays can be made 
before the CP asymmetries are measured. Regardless, it is clear that measurements of CP 
asymmetries and penguin decays will give complementary information. And in fact, unless 
the new particles are discovered in future colliders, it will be necessary to appeal to such 
measurements to infer their existence indirectly. 

The paper is organized as follows. In Section 2 we review and update the contributions 
of various models of new physics to B-B mixing, and hence to CP asymmetries in B decays. 
We summarize the current experimental constraints on the new-physics parameters which 
determine these contributions. For those models which can affect the SM predictions 
for the CP asymmetries, in Section 3 we examine their contributions to flavor-changing 
penguin decays. We conclude in Section 4. 

2. B-B Mixing and New Physics 

There are a variety of models of new physics which can contribute to B®-B® mixing 
(q = <i, s), and which therefore can affect CP asymmetries in B decays. In this section 
we review and update the contributions of these models to this mixing. (Note that we 
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include several models not discussed in Ref. [13].) We also examine how the new-physics 
parameters are constrained by current experimental data. In all cases, we search for new 
contributions to B®-B® mixing at least comparable to that of the SM [16]: 



*\2 



(3) 



where the mass difference AM is related to Mi 2 by AM q = 2|Mi 2 (-B ? )|, x t = m 2 /M, 
and 

3 1 3 x 2 \nx 



Mx) = 



1 9 1 

+ 



4 41 -x 2(l-x) 2 2(l-xY 



(4) 



New-physics contributions to B®-B® mixing are constrained by the measurements of 
the neutral S-meson mass differences, AMd [17] and AM S [18]: 



AM d = (0.470 ± 0.017) ps' 1 , AM S > 7.8 ps~ l 



(5) 



These values are consistent with the SM prediction [Eq. (3)], and constrain the CKM 
elements V t d and V ts as follows [2]: 



0.15 < 



V td 



< 0.34 



V, 



ts 



v cb 



> 0.6 



(6) 



These limits include the experimental errors on mt and V c b, as well as the theoretical error 
on f Bq \JB Bq . The bounds on these quantities due to the unitarity of the 3x3 CKM matrix 
alone are 



0.11 < 



V td 



v cb 



< 0.33 , 



V, 



ts 



v cb 



(7) 



With the addition of new contributions to B®-B® mixing, the constraints of Eq. (6) on 
Vtd and Vt s are relaxed, although the degree of relaxation is model-dependent. In certain 
models, the CKM matrix remains unitary, which implies that the bounds of Eq. (7) still 
hold. In other models, the 3x3 CKM matrix is not unitary, so that V t d or V ts can be 
much smaller and in principle even vanish, in which case B^-B^ or B®-B® mixing comes 
entirely from new physics. 

In order to see how large the new-physics contributions to B®-B® mixing can be in 
specific models, it is convenient to normalize these terms by the corresponding W box- 
diagram terms which appear in the SM, which are proportional to V 2 q . However, one should 
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note that in some cases the latter parameters can take values outside the SM constraints. 



by M-L2 rather than by Mf 2 M . 

2.1 ) Four generations [19] 

This is a model with an additional generation of quarks and leptons, including a new 
charge 2/3 quark, t' . The CKM matrix is 4 x 4, which can be parametrized by 6 angles 
and 3 phases. The unitarity triangle thus now becomes a quadrangle. There are new 
loop-level contributions, involving internal t' quarks, to both B®-B® mixing and penguin 
decays. The additional phases in the CKM matrix can play a role in the CP asymmetries. 

There is a model-independent lower bound of 45 GeV on the mass of the t' coming from 
LEP. There are stronger constraints on rrif of 0(100) GeV coming from hadron colliders, 
but these can be evaded since they depend on how strongly the t' couples to the b quark. 
There is an upper bound of 550 GeV on % coming from partial- wave unitarity [20]. A 
heavier t' will lead to a breakdown of perturbation theory. The strongest constraints on 
the CKM matrix elements involving the t' quark come from unitarity. There are additional 
constraints on the t' mass and its charged- current couplings coming from the K L -K S mass 
difference, from \e K |, from B^-B^ mixing, and from b — > sj. Since the measurements of all 
these observables agree with the predictions of the SM, they provide upper limits on their 
respective t' contributions, assuming no accidental cancellations. However, if one allows 
for such cancellations, the constraints become correspondingly weaker. Finally, we note 
that the fourth-generation neutrino must have a mass m v > M z /2 due to constraints from 
LEP. Since this is quite unlike the first 3 generations, many argue that the four-generation 
model is much less plausible. Still, it is a logical possibility. 

In this model, the extra phases in the 4x4 CKM matrix enter through the new 
contributions to B®-B® mixing. Assuming that this mixing is dominated by box diagrams 
with t and t' quarks, we have 



To avoid confusion in this respect, we will denote the W-hox contributions to B®-B® mixing 



M- 



-4— gen 
12 




G%M Bq n B Ml 
12n 2 



■f B B Bq [E(x t ,x t )(V tq V t t) 2 



(8) 



+ 2E(x t ,x t ,)(V tq Vf b )(V t , q V tlb .) + E{x t ,,x t ,){V t , q V t \) 2 ] , 



where 




(9) 
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Since the additional contributions to the mixing can be of a similar size to that of the SM, 
but with different phases, the CP asymmetries can be considerably altered. For example, 
the experimental value of B^-B^ mixing can be explained in the SM if V t d = 0.01 (for 
m t = 170 GeV, Vtb = 1); the phase of the mixing is then arg(VtdV t * b ) 2 . In a 4-generation 
model, in which the 3x3 CKM matrix is no longer unitary, this mixing can be dominated 
by the fourth generation: e.g. V t d ~ 0, V t 'd = 0.005, V tb = V t > b — 1/V%, and m t > = 480 
GeV. The phase of the mixing is then arg(Vt'dV^ b ) 2 , which may be quite different from 
the SM. 

2.2) Z -mediated flavor-changing neutral currents [21] 

In these models, one introduces an additional vector-singlet charge —1/3 quark, and 
allows it to mix with the ordinary down-type quarks. Since the weak isospin of the ex- 
otic quark is different from that of the ordinary quarks, flavor-changing neutral currents 
(FCNC's) involving the Z are induced. The Zbd and Zbs FCNC couplings, which affect 
B decays, are parametrized by independent parameters Udb and U sb , respectively, which 
contain new phases: 

A COS ts\y 

There are, however, constraints on the FCNC couplings coming from the process 
B — > n + n~X. The current experimental bound on the branching ratio of this process is 
[22] 

BR(B -> /j+ij-X) < 5 x 10" 5 , (11) 
while the contributions of ^-mediated FCNC's to this process are 

BR(B -> li+li-X) 



BR(B -> fj,i/X) 



|^d6| 2 + |^ sb | : 


I 


\v ub 


2 i p 
~ J - ps 


V c l 


>l 2 



(^) 2 + «f ITPw I 'J".t°; .2 ' ( 12 ) 



where g% = —1/2 + sin 2 ^^, g% = sin 2 # w , and F ps ~ 0.5 is a phase-space factor. The 
FCNC couplings U q b, q = d,s are then constrained to be 



U, 



qb 



< 0.044 , (13) 

Vcb 

or, taking \V cb \ = 0.0388 ± 0.0036 [2], 

\U qb \ < 0.0017 ±0.0002 . (14) 
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(Similar constraints can be obtained from the bound onB-> vvX [23].) 

The Z-mediated flavor-changing couplings U qb can contribute to B°-B G mixing: 



M? 2 (B q ) = 



V2G F M B m 



± fs q B Bq (u; b r . 



12 



(15) 



Recall that there can be new, independent phases in U qb . 

Comparing the contribution of this new physics to B° q -B^ mixing with that of the SM, 



we find 



v/2V 



AMY G F Mlx t f 2 (x t )\ Vtd V tb \ 2 
where we have taken \V t b\ = 1 and m t = 170 GeV. 
In this model the CKM matrix is not unitary: 



80 



\Vtd\ 



2 ' 



(16) 



v: d v ub + v; d v cb + v t * d v tb = u db , 



(17) 



so that the constraint of Eq. (7) on V td does not hold. Since \U db \ is bounded by Eq. (14), 
we find 

V td 



0.07 < 



Vcb 



< 0.37 . 



Consequently, 



AM? , 
——4- = 0.9-26 



\u db /v c 



cb 



0.04 



(18) 



(19) 



where the numerical coefficients 0.9 and 26 correspond to the largest and smallest values 
of iVtd/Kbl- The sum of the W and Z contributions to B d -B^ mixing is consistent with 
measurement for the entire range of Vtd- 

Bg-B® mixing can be analysed similarly. However, in this case, the effect of U sb on 
the violation of (sb) CKM unitarity is small, so that 



AM S Z 



0.15 



\u sb /v 



cb 



0.04 



(20) 



From this we see that B^-B^ mixing can in fact be dominated by Z-mediated FCNC. 
And although B° s -B® mixing is still mainly due to the W-hox contribution, the new-physics 
contribution may be non-negligible, so that the new phases in U sb can be important. 
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Thus, in both cases, measurements of CP asymmetries can differ considerably from the 
predictions of the SM. 

It is interesting to note that there exist specific models with seesaw-like predictions 
for the flavor-changing Z couplings [24]: 

^ = Vl|?> M = 0(0.1 — 1 TeV) . (21) 

Depending on the precise value of M, these couplings do not lie too far below the present 
limit [Eq. (14)], and may thus give sizeable contributions to B° d -B^ and B° s -B® mixing. 

In one particular flavor-changing Z model [25], CP is violated spontaneously, the 
CKM matrix is essentially real, and the new contributions to B®-B® mixing lead to very 
small phases. The unitarity triangle becomes a straight line and all CP asymmetries are 
expected to be tiny. 

2.3) Multi-Higgs-doublet models 

Models with more than one Higgs doublet can be classified into two types: (i) models 
with natural flavor conservation [26] , in which there are no flavor-changing neutral currents, 
and (ii) models in which flavor-changing interactions can be mediated by neutral scalars 
[27]. We discuss these in turn. 

(i) Natural Flavor Conservation 

In models with natural flavor conservation [28], the new charged scalars may give 
significant contributions to B®-B® mixing if their masses lie in the range of 50 GeV to 
about 1 TeV [29]. The Yukawa couplings of the charged scalars to up- and down- type 
quarks are given by 

C H ± = -J—J2 {^U L V CKM M D D R + Y l U R M u V CKM D L ) H+ + h.c. (22) 

Here U (D) is a vector of up- type (down-type) quarks, and M v (M D ) is the diagonal 
charge 2/3 (charge —1/3) quark mass matrix. Xj and V; are complex coupling constants 
arising from the mixing in the scalar sector. In the case of two Higgs doublets, in which 
one doublet provides masses to all quarks and the other decouples from the quark sector 
(model I), Y = —X = cot/?, where tan/3 is the ratio of the two vacuum expectation 
values. In a more popular version of the two-Higgs-doublet model (model II), found in 
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supersymmetric models for example (see discussion below), one scalar (0i) gives mass to 
the up- type quarks while the other scalar (</> 2 ) gives mass to the down- type quarks. In this 
case X = tan/3 and Y = cot (3 = v 2 /vi. 

For simplicity of presentation in the following we assume that only one charged Higgs 
is light; the others are heavy and decouple. This leaves two complex coupling constants, 
X and Y. There are several useful observations regarding Eq. (22). First, the Y term is 
dominated by the t-quark: C Y H± ~ Y (m t /M w ) H + i R V t id iL . Second, due to the smallness 
of the down-type quark masses, the X term is important only if |X| 3> \Y\. In this region 
of parameter space, CP violation can appear in charged-Higgs exchange if X and Y have 
a nonzero relative phase. However, the observed rate of b — > S7 constrains Im(XY*) < 2-4 
[30] , thus ruling out the possibility that CP-violating effects due to charged-Higgs exchange 
can compete with those due to W exchange. Since the inclusion of the X term does not 
lead to new CP violation, and since in general it is much smaller than the Y term, from 
here on we will generally ignore the X term altogether. We will refer to the possibility of 
very large X only when its effect is particularly important. 

There are two types of box diagrams involving charged Higgs bosons which contribute 
to -Bg-Sg mixing: those with one H and one W, and those with two _£Ts. The total 
charged-Higgs contribution is given by [31] 

Mf 2 + (B 9 ) = — ' 4 ffi w f B B Bq (V tq V t tr [I HH + I HW ] , (23) 



where 



I HH =x t y t I 1 (y t )\Y\ 4 , I HW = 2x t y t [AI 2 (x u y t )+h(x u y t )]\Y\ 2 , (24) 



with 

^) = (I^+(T^)3 . (25) 

.-r , x t / \ (x — 4-y)lnu 3x\nx x — A , nnS 

4I 2 (x, y) + I 3 (x, y) = / ^ y + ^ + — . 26 

{y-x)(l-y) 2 {y-x)(l-x) 2 (1 - x)(l - y) 

Here, x q = m^/M^ and y q = m 2 q /M^ + . Note that M" 2 (B q ) involves the same CKM 
factors and has the same phase as the VF-box contributions. Therefore, the ratio of the 
two terms is independent of CKM factors and is positive, such that the two contributions 
add up constructively. For M H + in the range 100-400 GeV, this ratio may be approximated 
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to within about 20% by an inversely linear relation: 



For smaller and larger Higgs masses this simplified expression holds within about 30%. 

The parameters \Y\ and M H + are constrained by comparing the observed rate of 
b — > S7 [32] with the SM prediction [33] . These constraints also depend on X [34] . At 3cr 
the bounds are 



-0.56 < \Y\hs w (y t ) + XY*G H (y t ) < 0.27 , (28) 



where 



G w {y) = 12(i y _^ )4 [(7 - 5y - 8y 2 )(l - y) + 6y(2 - Sy)ln(y)] 
G H {y) = 6( ^ 1)3 [(3 - 5y)(l -y) + 2(2 - 3y)ln(y)] . 



(29) 



The implications of these bounds on Y and M H + depend somewhat on the details of 
the model [35]. In a two-Higgs-doublet model of type II (XY = 1), charged-Higgs masses 
below about 300 GeV are already excluded, independent of the value of Y . In a two-Higgs- 
doublet model of type I (X = —Y), Higgs masses in the entire range M H + < 800 GeV are 
excluded if one assumes \Y\ > 2.7. However, in a general multi-Higgs-doublet model, in 
which X and Y are independent parameters, the constraints become weaker. This leaves 
a large region of \Y\-M H + parameter space in which the charged-Higgs contribution to 
Bq-B® mixing [Eq. (27)] can be significant and even dominant. For example, in a general 
model, the values M H + = 400 GeV and Y = 3 are allowed, which implies Mf 2 + /M^ = 7. 
In this case the sum of the W and H + terms (dominated by H + ) is consistent with the 
measurement of B®-B® mixing for the smallest values of Vtd in the unitarity range [Eq. (7)]. 

There is also a bound on the parameters \Y\ and M H + from the latest ALEPH mea- 
surement of R b = F(Z -> bb)/T(Z -> hadrons) = 0.2158 ± 0.0014 [36]. Assuming the 
neutral-Higgs contribution to Rb is small (i.e. X is not too large) [37], at 3a this results 
in the constraint [38] 



where y t = m^/M^ + and 



\Y\ 2 F(y t ) < 1.7 , (30) 



F ^) = (^1)2 fo-l-lny]. (31) 



This constraint is somewhat weaker than that from b — > s"f. 
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It is important to note that the phase of B®-B® mixing is unaffected by these new 
contributions. Also, the unitarity of the 3x3 CKM matrix holds in these models. Con- 
sequently, although the extraction of \V t d\ through the measurement of AMj has to take 
the H + contribution into account, the measurement of its phase in the asymmetry of 
B d — > ^K s is unaffected by the presence of the new physics. 

Particularly interesting are models of spontaneous CP violation, in which the entire 
Lagrangian is CP invariant while the vacuum is not. (In the Weinberg three-Higgs-doublet 
model [26], this possibility seems to have already been ruled out by the experimental 
upper limit on the neutron electric dipole moment [30].) In this case CP is violated in 
(neutral and charged) Higgs exchange, while natural flavor conservation leads to a real 
CKM matrix [39]. Thus, the unitarity triangle becomes a straight line, and the amplitude 
of B®-B® mixing is real. 

(ii) Flavor-Changing Neutral Scalars 

In the second class of models flavor- changing neutral scalar interactions between 
quarks i and j exist, but are suppressed by factors Fij due to an approximate global 
symmetry. Such models have recently received special attention [40]- [44]. In this case 
B®-B® mixing may also receive large contributions from tree-level neutral Higgs exchange 
amplitudes which carry new phases. Denoting the flavor-changing neutral Higgs couplings 
by (rrii/v)Fij (m.j > m,j, v~ 2 = \^2G F ), their contributions to B®-B® mixing are given in 
the vacuum insertion approximation by 

M 12 (B q ) « — t qh , (62) 

where M H o includes possible complex mixing among several neutral Higgs fields. 
Comparing with mixing in the SM, we find 

^H^^V (™GeVV (33) 



Thus the neutral Higgs contributions may substantially modify the SM prediction for B®- 
B® mixing. The unitarity triangle holds in this model. However neither the magnitude of 
V t d nor its phase can be directly measured through AM d and the asymmetry inBj — > ^K s , 
respectively. 

In models with specific predictions for F^, this leads to large effects in an interesting 
range of Higgs masses. For instance, neutral-Higgs contributions to B®-B® mixing are 
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sizeable for Higgs masses around 100 GeV when flavor-changing couplings have the form 



and for Higgs masses of a few hundred GeV up to about a TeV in models in which [40] 



In general, these new contributions carry unknown phases and have to be added to the 
charged Higgs contributions [Eq. (27)]. 

In the presence of flavor-changing scalar interactions, the special case of spontaneous 
CP violation does not, in general, forbid a phase in the CKM matrix. For a particular 
choice of the softly-broken symmetry this phase may, however, be very small [44] or may 
even vanish [42]. This would imply that the unitarity triangle becomes a straight line, 
while the B®-B® mixing amplitude carries a complex phase. 

2.4) Left-right symmetric models 

In left-right symmetric models [45], the gauge group is extended to SU(2) L x SU(2) R x 
U(1) B _ L: along with a discrete L <-> R symmetry. The right-handed CKM matrix is then 
related to its left-handed counterpart: V R = V L or V R = (V L )*. The right-handed W R 
can participate in weak processes in the same way as the ordinary W (although certain 
decays are forbidden if the v R is too heavy). In particular, the W R can contribute to 
Bq-B® mixing through box diagrams, as in the SM. However, limits from the K L -K S mass 
difference constrain the W R to be heavier than 1.4 TeV [46], which would render its effects 
in the B system negligible. 

If one abandons the discrete L <-> R symmetry, so that V R is unrelated to V L , the 
constraints from the K L -K S mass difference can be evaded. For judicious choices of the 
form of V R , CP-conserving experimental data permit the W R to be considerably lighter, 
M R > 300 GeV [47]. However, unless the elements of V R are considerably fine-tuned, there 
will be large contributions to the CP- violating parameter \e K \ [48]. Assuming no such fine 
tuning, the W R is again constrained to be heavy, M R > 5 TeV. 

The above analysis assumes that V L has the same form as in the SM. However, this 
need not be the case. For example, it was suggested [49] that B decays might in fact 
be mediated by the W R , instead of the ordinary W. The long B lifetime would then 
be interpreted as being due to the heaviness of the W R: rather than to the smallness of 



[41], [42] 



F qb = V tq V t 



(34) 




(35) 
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V c t,. A variety of different forms for V L and V R were proposed [49], [50]. In all cases, 
the V* d element was considerably smaller than in the SM, leading to the prediction that 
BR(b — > ccd)/BR(b — > ccs) is at most O(10~ 4 ). However, the decay B — > *7r has since 
been observed with a branching ratio in agreement with the SM [51], effectively ruling out 
all such models. 

Our conclusion is therefore that there are no important new-physics effects in the B 
system within left-right symmetric models. The one possible exception is if one consider- 
ably fine-tunes the right-handed CKM matrix [52] , but we do not consider such possibilities 
here. 

2.5) Supersymmetry 

In the supersymmetric standard model (SSM) [53], the gauge group is unchanged, 
but a plethora of new particles is added. These include the supersymmetric partners of 
the SM particles, as well as a second Higgs doublet (in some versions, additional Higgs 
representations are also present) . In the SSM there are thus a variety of new contributions 
to B®-B® mixing. These come from box diagrams with internal (i) charged Higgs bosons 
and charge 2/3 quarks, (ii) charginos and charge 2/3 squarks, (iii) gluinos and charge —1/3 
squarks, and (iv) neutralinos and charge —1/3 squarks. The relative sizes of these new 
contributions, as well as their phase information, depend on the version of the SSM. 

We first consider the minimal supersymmetric standard model (MSSM), which is 
usually taken to mean the low-energy limit of the minimal spontaneously-broken N = 1 
supergravity model. Here one typically imagines that there is unification at some high scale 
(M x ), and that supersymmetry (SUSY) is broken at this scale by some unknown mecha- 
nism (e.g. a "hidden sector") which interacts only gravitationally with the known fields. 
SUSY breaking is parametrized by soft breaking terms in the supergravity lagrangian [54] . 
The low-energy effects of SUSY breaking, as well as the masses of the superpartners, are 
calculated by running the renormalization group equations down from M x to the weak 
scale. The net effect is that, in addition to the usual gauge and Yukawa couplings, the 
MSSM is described by only 4 new parameters. The masses and mixings of all superpart- 
ners at low energy can be described in terms of these 4 parameters. If one also requires 
that the spontaneous breaking of SU(2) L x U(1) Y be induced radiatively, there is a further 
reduction in the number of SUSY parameters from 4 to 3 [55] . 

Although there are several new SUSY contributions to B®-B® mixing, in the MSSM 
these all have the same phase as in the SM, to a good approximation. We consider them 
in turn [56]: 



16 



• (i) charged Higgs bosons and charge 2/3 quarks: These contributions have already been 
described above in the discussion of the two-Higgs-doublet model (Sec. 2.3 (i)). Unless one 
goes to extremely large values of tan/?, the charged Higgs couples the down-type quarks 
only to the t-quark: 

£ H ± ~ 4=cot/? ^H + i R V ti d iL . (36) 

Thus the contribution to B®-B® mixing from charged Higgs bosons is proportional to 

(V t * b V tq ) 2 , as in the SM. 

• (ii) charginos and charge 2/3 squarks: The couplings of the and to up-type 
quarks are very similar to those of the and H^. (Note that the physical charginos 
are in general linear combinations of and -fi .) In particular, the couples only to 
left-handed up-type squarks: 



C~ =4=(u, 5, t) L V CKM W lL hi, (37) 



V2 




while in the limit of negligible down-type quark masses, the couples mainly to right- 
handed squarks (assuming non-extreme values of tan j3) : 

C ^f2^W hV " llLd " (38) 

where 7^ = (1 — 7s)/2. The contributions to B®-B® mixing of both the W ± and H ± are 
proportional to (V t * b V tq ) 2 , as in the SM. For the H ± this follows directly from the above 
equation, while in the case of the W ± , one uses the unitarity of the CKM matrix, along 
with the fact that rriu L = rric L 7^ in the MSSM, to arrive at this result. 

• (Hi) gluinos and charge —1/3 squarks: The important coupling of the gluino (g) to 
down-type quarks and squarks is 



Cg = V2g 3 (d, s, b) L VcK M giL \ s . (39) 




(There is also a generation-diagonal coupling involving right-handed down squarks, but 
this cannot contribute to B^-B® mixing.) Note that the coupling is proportional to V CKM . 
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This, along with the fact that mj^ = m§ L ^ in the MSSM, leads to a contribution to 
B®-B® mixing proportional to (V t * b V tq ) 2 , as in the SM. 

• (iv) neutralinos and charge —1/3 squarks: The physical neutralinos are linear combina- 
tions of the photino (7), the Zino (Z), and the two neutral Higgsinos (#12). The couplings 
of the 7 and Z to down-type quarks and squarks are similar to that of the g [Eq. (39)]: 



d 

e ) (d, s, b) L V CK MllL I s 



£~ = - 

z cos 9,, 



-- + - sin 2 6L 
2 3 



d 

(d, s, b) L V CKM Z^ L I s 



(40) 



The dependence on V CKM of the couplings of the 7 and Z is just like that of the gluino, 
leading to a contribution to B®-B® mixing which is proportional to (V t * b Vt q ) 2 . As for 
the neutral Higgsinos, their coupling to down-type quarks and squarks is proportional to 
M D /M W cos/3, which is negligible (unless one goes to extremely large values of tan/3). 

Thus, to a good approximation, in the MSSM all new SUSY contributions to B^-B^ mixing 
have the same phase as in the SM. Therefore the CP asymmetries in B decays will not 
be modified. The full expressions for these contributions are quite complicated, so we do 
not reproduce them here (we refer the reader to Ref. [56]). The strongest constraints on 
the SUSY parameters come from direct searches. For example, there are lower bounds of 
176 GeV and 45 GeV on squark and chargino masses, respectively [57]. The parameter 
space is sufficiently complicated that it is very difficult to establish firm constraints from 
loop-level processes such as b — > S7. We note, however, that the effect of supersymmetry 
on -B°--Bg can be quite significant - for certain values of the parameters, the total SUSY 
contribution to B%-B% mixing can be twice as large as that of the SM. 

Recently [58] it was suggested that CP violation could possibly come from SUSY 
breaking alone, with a real CKM matrix. In this case, which is essentially a superweak- 
type model, the unitarity triangle becomes a straight line and all CP asymmetries in B 
decays are very small. 

We now turn to nonminimal SUSY models. For generic squark masses, supersym- 
metric contributions enhance flavor-changing processes such as K°-K° well beyond their 
experimental values. Any nonminimal SUSY model should address this problem (in the 
MSSM this problem is resolved since, to a good approximation, the squarks are degener- 
ate). 
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One class of nonminimal SUSY models solves the FCNC problem by imposing an 
Abelian horizontal symmetry on the lagrangian [59] . In this case the quark mass matrices 
are approximately aligned with the squark mass-squared matrices. This has the effect 
that the mixing matrix for quark-squark-gluino couplings is close to a unit matrix, so that 
FCNC are suppressed, even though the squarks may not be degenerate. In most such 
models, the SUSY contributions to B®-B® mixing are quite small. However, it is possible 
to construct models in which Mf 2 ysy (B^)/Mf 2 M (.E?2) is as large as 0.15, with a negligible 
effect on B®-B® mixing. Since the phase of the SUSY contribution is unknown, this can 
lead to measurable deviations from the SM predictions in CP asymmetries involving B^ 
decays. In another class of nonminimal SUSY models, known as effective supersymmetry 
[60], the suppression of FCNC's applies only to the first two families of squarks. In this 
case new-physics effects on B®-B® mixing can be much larger [61]. 

Another approach which is often taken is to ignore the FCNC problem altogether. One 
simply assumes that all SUSY parameters take the maximum allowed values permitted by 
experiment. In such "models" one can have non- negligible contributions to B®-B® mixing 
which have different phases than in the SM. These contributions typically involve right- 
handed squarks. For example [62] , the general coupling of a H ± to down-type quarks and 
up-type squarks can be written 



C H^^ MwSin ^ ~ c ' t)nU^U u R M u V CKM H lL [ s ] , (41) 




where U% (U%) is the transformation matrix of the right-handed up-type squarks (quarks) 
needed to diagonalize the squark (quark) mass matrix. In the MSSM, U% = U%, leading to 
a contribution to B®-B® mixing proportional to (V* b Vt q ) 2 [Eq. (38)]. However, in general 
this relation need not hold, in which case there can be new phases in this contribution to 
By-B® mixing. 

As another example, consider again the contribution of gluinos and charge —1/3 
squarks to mixing. In the MSSM there is no intergenerational mixing among 

right-handed down-type squarks. As a consequence, the contribution to B^-B^ mixing of 
gluinos and charge —1/3 squarks involves only left-handed squarks [Eq. (39)]. However, in 
nonminimal SUSY models, this need not be the case [63] - there can be intergenerational 
mixing among right-handed down-type squarks. In general, this mixing matrix is unrelated 
to V C kmi so that there can be new phases in this contribution to B^-B^ mixing. 

The main problem with this approach is that there is little predictivity. There are, 
in general, a very large number of parameters - the masses of the superpartners, their 
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mixings, etc. Thus, although one can describe how new phases can enter B®-B® mixing, 
it is virtually impossible to analyse such effects in a systematic way. 

3. Penguin Decays and New Physics 

As discussed in the introduction, any new physics which contributes to B®-B® mixing 
will also contribute to flavor-changing B decays. Before examining the new-physics contri- 
butions to such penguin decays, we first review the SM predictions. Two aspects of these 
predictions are of particular interest to us: (i) the actual size of the branching ratios for 
various penguin decays, and (ii) the uncertainties, both experimental and theoretical, of 
the predictions. New-physics effects will be considered important in a particular penguin 
decay only if they change the branching ratio by quite a bit more than the uncertainty 
in the SM prediction - in other words, we are looking for "smoking gun" signals of new 
physics in such decays. 

3.1) The standard model 

• b — > qj, q = d,s: The lowest-order amplitude for the decay b — > q^ is [33] 

A(b - qi ) = nVidFiix^e'sv^im^ + m slL )b , (42) 

v i 

where the sum is over the up-type quarks, and q^ and are the photon's four-momentum 
and polarization, respectively. The function F2 is given by 

F 2 (x) = X _ [6x(3x - 2) logo: - (x - l)(8x 2 + 5x - 7)] . (43) 

Due to the smallness of the u- and c-quark masses, F2(x u ), F2(x c ) <^ F2(xt), so that the 
b — > (77 amplitude is dominated by t-quark exchange. A full quantitative treatment of 
these decays requires the calculation of the important QCD corrections. Including these, 
the SM branching ratios are [33] 

BR(B -> X s 7) = (3.2 ± 0.58) x 10" 4 , 

' (44) 

BR{B ^X dl ) = (1.0 ± 0.8) x 10" 5 . 

The uncertainties include both experimental errors [m t , B semileptonic branching ratio] 
and theoretical errors [fx (the renormalization scale), A QCC , and the ambiguity in the 
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interpretation of m t (pole or running mass)], combined in quadrature. For b — > d'j, 
\V t d/V ts \ = 0.24 ± 0.11 has been used and combined in quadrature with the other errors. 
The decay b — > S7 has actually been measured by the CLEO collaboration [32] : 

BR(B -> X s 7) = (2.32 ± 0.67) x 10" 4 . (45) 

This measurement can be used to constrain models of new physics, as already demonstrated 
in Eq. (28). 

• b — > ql + l~ , q = d, s: This class of decays is rather complicated theoretically [33]. First, 
one must calculate, at next-to-leading order, the Wilson coefficients of ten local operators 
which mix under renormalization. Second, one needs the matrix elements relevant for 
B — > X S) dl + l~, which can be calculated using the spectator model, along with (3(l/m 2 ) 
corrections. Finally, long-distance effects due to J/ty and resonances must also be taken 
into account. The short-distance contribution for b — > sl + l~, which can be measured far 
away from the resonances, gives the following branching ratios, taken from Ref. [33]: 

BR(B -> X s e+e") = (8.4 ± 2.2) x 10" 6 , 

BR(B -> X s //+//") = (5.7 ± 1.3) x 10" 6 , (46) 
BR(B -> X s r+r") = (2.6 ± 0.5) x 10" 7 . 

For b — > dl + l~ , the branching ratios are 

BR(B -> X d e+e") = (4.9 ± 4.3) x 10" 7 , 

BR(B -»• X d n + n~) = (3.3 ± 2.8) x 10" 7 , (47) 
BR(B -> X d r+r-) = (1.5 ± 1.3) x 10" 8 . 

For all decays, the errors come from the same sources as in b — > 57: mt, B semileptonic 
branching ratio, the renormalization scale /x, A QCD , and |Vfd/Vt s | = 0.24 ±0.11. 

• B® — > / + /~, g = d, s: This decay can be calculated quite precisely in the SM. By 
including the QCD corrections [64], the renormalization-scale uncertainty is reduced to 
(9(1%). There is still some hadronic uncertainty, parametrized by the S-meson decay 
constant f B . The branching ratios are [33] 

BR(B° S -> r+r") = (7.4 ± 2.1) x 10" 7 (/ Ss /232 MeV) 2 , 
BR(B° S -> = (3.5 ± 1.0) x 10" 9 (/ Ss /232 MeV) 2 , 

BR(B° d -> r+r") = (3.1 ± 2.9) x 10" 8 (/ Sd /200 MeV) 2 , 



BR(B G d -> ^+^") = (1.5 ± 1.4) x 10" 10 (/ Sd /200 MeV) 



2 
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(The branching ratios to e + e~ are some 5 orders of magnitude smaller than those for 
.) The error in the B® branching ratios is due to the uncertainty, both experimental 
and theoretical, in the top-quark mass. The B% branching ratios have a larger error due 
to the Vtd CRM matrix element: |Vfd/Vts| = 0.24 ±0.11. At present, the best upper limits 
are BR(B° S -> n + n~) < 8.4 x 10" 6 and BR(B° d -> n + n~) < 1.6 x 10" 6 [65], with no 
significant limits on the t + t~ final state. 

• Gluon-mediated exclusive hadronic decays: These arise from the quark-level process b — > 
qq'q'. Throughout this paper we will refer to such loop-level decays as "hadronic penguins." 
There are two ingredients needed to calculate the rates for hadronic penguin decays in the 
SM. First, the rates for the quark-level decays b — > sqq and b — * dqq are computed. This is 
done similarly to the decays b — > s and b — > dl + l~: the Wilson coefficients of a variety 
of operators are calculated as one renormalizes down from the weak scale to the b- mass [66] . 
This can be done with reasonable precision. Second, one calculates the hadronic matrix 
elements for the hadronization of the final-state quarks into particular final states [67]. 
It is this step which introduces enormous uncertainty. These hadronic matrix elements 
are typically evaluated using the factorization approximation. Unfortunately, it is difficult 
to estimate the error incurred by applying this approximation to penguin decays. The 
predicted rates for exclusive hadronic penguin decays can easily be in error by a factor of 
2 to 3. (Much of this uncertainty cancels in the ratio of rates of corresponding b — > d and 
b — > s processes, which is given in the SM by IV^/V^I 2 [68].) 

Since the SM predictions for hadronic penguins have considerable uncertainties, if one 
wants an unmistakable signal of physics beyond the SM in such decays, the new-physics 
effects must be enormous - they must change the SM rates by an order of magnitude or 
more. It is therefore sufficient for our purposes to obtain approximate, order-of-magnitude 
estimates for both the SM and new-physics effects. To this end, we will use the following 
approximate form for the amplitude of the SM gluonic penguin contribution to the decay 
b — > qq'q': 

A penguin ~ log (^j V tq V t t ~ 0.04 V tq V t l . (49) 

(Note that the coefficient 0.04 is about the same size as the largest of the Wilson coefficients 
of penguin operators [66].) 

This expression can be used to estimate the order-of-magnitude rates for b — > s and 
b — > d penguins in the SM. For example, the branching ratio for B d — > which is 

dominated by the tree-level b — > uud amplitude, is O(10 -5 ). Comparing this decay with 



22 



b — > s penguins, which dominate B^ — > 7r we find 



-^penguin 






®MV t lV ts 




uud) 




v* b v ud 



where we have used \V ts \ ~ |V c b| and iV^b/Kbl = 0.08. This is consistent with the ob- 
servation of a combined sample of B® — > 7r + 7r~ and — > 7r~K + decays [69], in which 
about equal mixtures of both modes are most likely. Thus, assuming that the hadronic 
matrix elements of tree and penguin operators have similar magnitudes, we expect pure 
b — > s penguin hadronic decays (e.g. B + — > n + K s , B^ — > (pK s ) to have branching ratios 
of O(10 -5 ). b — > d penguins can be analyzed similarly: 



Ap en guin(b > C?) 




o.o4v;* 6 y td 









(51) 



Since 1.4 < |Vtd/Kjb| < 4.6 [2], the penguin amplitude is about 1/10 the size of the 
tree amplitude, and may be larger if the hadronic penguin matrix elements are enhanced 
relative to those of tree amplitudes. Thus, pure b — > d penguin hadronic decays (e.g. 
B + — > K + K s ) should have branching ratios of O(10 -7 ), or somewhat larger. 

As an aside, we note that Eq. (51) demonstrates why penguin contamination is a 
concern in the extraction of sin2a using the CP asymmetry in B^ — > n + n~ . According to 
this estimate, the penguin amplitude can be as much as ~ 15% of the tree amplitude in 
magnitude (or even larger, if penguin matrix elements are enhanced), and with a different 
phase. This can lead to considerable uncertainty in the extraction of sin2ct, and shows 
why isospin techniques [5] are necessary to remove the penguin contamination. 

• Electroweak-penguin- dominated exclusive hadronic decays: Here, the diagrams contribut- 
ing to the quark-level process b — > qq'q' consist of 7 and Z penguins and box diagrams. 
Of these, the diagram involving Z exchange is the most important since it is enhanced 
by a factor of m^/M^. Throughout this paper we will refer to such processes as "elec- 
troweak penguins." As with hadronic penguins, the calculation of SM rates for exclusive 
electroweak penguin decays suffers from large uncertainties in the hadronic matrix ele- 
ments. We will therefore again use an approximate form for the amplitude of the SM 
electroweak penguin contribution to the decay b — > qq'q'. This can be obtained from 
Eq. (49) through the replacements a s (mj) — * a^C^b) and log (mf/m^) — > mf/M^,. There 
is an additional factor of 2 due to a larger Wilson coefficient for the electroweak penguin 
operator. Therefore, 

A EWP ~ V tq V t t ~ 0.008 V tq V t t . (52) 
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Note that, if one combines the final-state q and q' quarks to form a meson, there is an 
additional color-suppression factor, a 2 [70]. Compared to color-allowed decays (i.e. forming 
a meson from q' and q'), which are parametrized by ai, this suppression is a^/ai = 0.2 
[71]. 

Comparing Eqs. (49) and (52), we note that electroweak penguin amplitudes are 
suppressed relative to their hadronic penguin counterparts by a factor 0.2. Therefore we 
expect b — > s and b — > d electroweak penguin decays to have branching ratios in the range 
10~ 7 -10~ 6 and 10~ 9 -10~ 8 , respectively [72]. Some examples of decays which are dominated 
by electroweak penguins are B® — > (pn (b — > s) and B + — > (pn + (b — > d). 

At this point several observations are in order. From the above summary, we see that 
the SM b — > d penguins have branching ratios which are about one to two orders of 
magnitudes smaller than their b — > s counterparts. Therefore, unless the new physics has 
a large influence on the b — > d FCNC, its effects on b — * s penguins will be detected first. 
On the other hand, CP asymmetries involving decays are likely to be measured — and 
hence will reveal the presence of new physics — well before those involving B® mesons. So, 
from a practical point of view, this poses a bit of a problem. That is, even if new physics 
is detected in CP violation in B^ decays, it will be possible to test its nature by looking 
at b — > d penguin decays only if the small SM branching ratios for these processes are 
significantly enhanced. Conversely, if one finds new-physics effects in b — + s penguins, it 
will be difficult to determine its origin by looking at CP asymmetries in B® decays. 

This having been said, however, the situation is not quite so bleak. Although it is 
possible to construct models of physics beyond the SM in which only one of the b — > d 
or b — > s FCNC's is changed, in practice both FCNC's are affected in most new-physics 
models. For example, in models with Z-mediated FCNC's, the b — ► d or b — »• s FCNC's are 
indeed described by different parameters. However, both of these FCNC's arise due to the 
mixing of the charge —1/3 quarks with an exotic vector singlet quark. It would be difficult 
to imagine that only one of the two FCNC's is induced, although this is a logical possibility. 
As another example, consider models with four quark generations. The enlarging of the 
CKM matrix to 4 x 4 will, in general affect both FCNC's. Thus, in looking for new physics, 
the measurements of CP asymmetries and B penguins are complementary, and it is quite 
likely that, should new physics be discovered, its nature will be revealed only by studying 
both B asymmetries and decays. 

We now turn to the new-physics contributions to the above penguin decays. 
3.2) Four generations 
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Models with four generations have a number of new parameters which can enter in 
B decays: m t > and the CKM matrix elements involving the t' quark. Furthermore, the 
strongest constraints on V ts and V t d in the SM come from the unitarity of the 3x3 CKM 
matrix. When this matrix is enlarged to be 4 x 4, these constraints are weakened, so that, 
in effect, V t d and V ts are also unknown parameters. 

The parameter space of four-generation models is therefore quite large. Rather than 
exploring the entire space, we will simply present an "existence proof." That is, we will 
show that it is possible to choose values of m t ' and the CKM matrix elements, consistent 
with experimental data, which significantly affect both CP asymmetries and B penguin 

decays. 

As discussed in Sec. 2.1, the experimental value of B^-B^ mixing can be reproduced 
if V td ~ 0, V t 'd = 0.005, V tb = Vt> b ~ l/y/2, and m t < = 480 GeV. The phase of this 
mixing may, however, be quite different than in the SM. In this scenario, all penguin 
decays involving the b-d FCNC will also be dominated by the fourth generation. Below we 
consider the effects of this choice of parameters on b — > d penguins. 

We first consider the decays b — > q^. The experimental measurement of b — > S7 can be 
easily accommodated by adjusting the t and t' contributions. However, the decay b — > dj 
will involve only t' exchange. Using Eq. (42), we find that \V t 'dV t 1 b F 2 {x t ') \ = 0.001, as 
compared to the SM value of \V t dV t * b F 2 (x t ) \ = 0.002. Therefore, this choice of parameters 
will result in a branching ratio for b — > dj which is roughly 4 times smaller than in the 
SM. Given the large uncertainties in the SM prediction, this cannot be considered an 
unmistakable signal of new physics. 

Now consider hadronic penguins. From Eq. (49), we have 
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(53) 



Thus, in this model, the branching ratios for exclusive b — > d hadronic penguins will be a 
factor of ~ 5 smaller than in the SM. Given the uncertainties in the SM predictions this 
is also a marginal "smoking gun" signal of new physics. 

The situation is better for electroweak penguins. From Eq. (52), 
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The branching ratios for exclusive b — > d electroweak penguins are thus a factor of 



(54) 
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larger than in the SM. This same enhancement applies to the decays — > l + l . In both 
cases, this would be a quite convincing signal of physics beyond the SM. 

Thus, in models with four generations, we have shown that there are regions of pa- 
rameter space in which both B CP asymmetries and B penguin decays are significantly 
affected. If a discrepancy with the SM is found in the measurement of the CP asymmetries, 
the study of the decays can help pin down the new-physics parameters. Admittedly, in the 
particular example we have chosen, the branching ratios for the affected processes are all 
small, of O(10~ 7 ) or smaller. But the key point here is that the various penguin processes 
depend differently on the masses of the t and t' quarks. It is thus straightforward to find 
a set of parameters, consistent with current experimental data, in which the b — > s FCNC 
decays are affected. In this case it would be CP asymmetries involving B° s decays which 
would be altered. Of course, since the four-generation CKM matrix is 4 x 4, in the general 
case both the b — > d and b — > s FCNC's will be changed from the SM, affecting all CP 
asymmetries and penguin decays. 

3.3) Z-mediated FCNC's 

For the process b — > qff, the amplitude due to Z-mediated FCNC's is 

AG 

M z - F cnc =^=-U qb q^ L bf [glnftfL + glltf*] f ■ (55) 
The rate for b — > qff is simply 

nb^«ff)=^\U, b \°\(d) 2 +( g i) 2 }, (56) 



while the rate for B® — > l + l is 



r(S; - l+r) = ^T Bq f B M Bq ml\U qb \ 2 . (57) 

With these expressions in hand, we can now calculate or estimate the contributions to 
penguin decays due to Z-mediated FCNC's. In all cases, when presenting numbers, we use 
the upper limit U qb < 0.0017 [Eq. (14)]. 

The upper limit of U qb < 0.0017 is in fact derived from the experimental limit on the 
branching ratio of B — > /i + \i~X [see Sec. 2.2]. Thus, with this value of U qb , Z-mediated 
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FCNC models "predict" that BR(B -> ^i+^X) = 5 x 10" 5 . For b -> s FCNC's, this is 
roughly an order of magnitude larger than the SM prediction, while for b — > d transitions it 
is about 2 orders of magnitude larger. If the branching ratios for the decays B — > X s l + l~ 
and £? — > Xal + l~ are observed to be consistent with the SM predictions, this will imply 
that \U s b\ < 6 x 10~ 4 and \Udb\ ^ 1 x 10 -4 . In both cases, this implies that the new- 
physics effects in B®-B® mixing are negligible. Conversely, if B®-B® mixing is significantly 
affected by ^-mediated FCNC's, then one expects to see a substantial enhancement of the 
branching ratios for B — > X s l + l~ and/or B — > Xd l + l~ . 

We now turn to B® — > Z + /~, q = d, s. From Eq. (57), the rates for these processes, 
due only to ^-mediated FCNC's, are 

BR(B° S -> t + t-)\ z ^ fcnc < 1.6 x 10" 5 (/ Bs /232 MeV) 2 , 
BR(B° S ^fi+fi-)\ z _ FCNC < 5.8 x 1Q~ 8 (f Bs /232 MeV) 2 , 

(58) 

- t+t-)| z _ fcjvc < 1.2 x 10" 5 (/ Bd /200 MeV) 2 , 
B#(B° -> ^/OU-^c < 4.2 x 10" 8 (/ Bd /200 MeV) 2 . 



Thus, if the FCNC parameters have values near the present upper limits, the predicted 
rates for B° s — > l + l~ and B® — > Z + /~ are respectively about 20 and 300-400 times larger 
than those expected in the SM. 

Turning to hadronic and electroweak penguins, we note that there are three types of 
comparisons which can be made. First consider decays such as B® — > K S K S (b — > s) or 
— > K S K S (b — > d), which receive contributions from ordinary (gluonic) penguins and 
color-suppressed Z-mediated FCNC's. For these decays we have 
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(59) 



where we have taken |f7 g b/V C 6| < 0.044 and | Vtd/V^b| ~ 0.24 in our estimates of the ratio. 
The branching ratios for this type of B decays will therefore not be significantly affected 
by Z-mediated FCNC's. 

Next we have decays such as B% — > <f>K s (b — > s) or B® — > (pK s (b — > d). Here there 
are contributions from ordinary gluonic penguins and color-allowed Z-mediated FCNC's. 
Then 
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(60) 
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In this case, Z-mediated FCNC's will not much affect the b — > s penguin decays, but 
the branching ratios for b — > d penguins can be increased by a factor of 3-4. Given the 
uncertainties in the SM predictions for such decays, this cannot be considered significant. 
However, if the value of V td is in fact smaller than 0.24 (in this model it can be as small as 
0.07), then the branching ratios for b — > d penguins will be correspondingly increased. We 
consider this to be a marginal prediction, since it depends sensitively on the true value of 

v td . 

Finally, decays such as — > (fin (b — > s) or B + — > (fin + (b — > d) have contributions 
from ordinary electroweak penguins (but not gluonic penguins) and Z- mediated FCNC's. 
Here 
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(61) 



The effects of Z-mediated FCNC's on such decays are clearly enormous. The branching 
ratios for pure electroweak penguin decays can be increased by as much as a factor of ~ 25 
(b — > s) or ~ 500 (6 — > d)l Clearly this is a "smoking gun" signal of new physics. 

On the topic of hadronic penguin decays, there is another possibility which should 
be mentioned. In the SM, the decay B + — > 7r + 7r° occurs principally via a tree-level 
amplitude - there is no gluonic penguin and the electroweak penguin is much suppressed. 
One therefore does not expect to find CP violation in this mode. However, Z-mediated 
FCNC's also contribute to this decay. Comparing this new-physics contribution with the 
SM tree-level amplitude, we find 
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< 0.5 



(62) 



The Z-mediated FCNC contribution to this decay could therefore be substantial. If there 
is a significant strong phase difference between the two amplitudes, following for instance 
from different rescattering in the 1 = 2 channel, this could lead to direct CP violation in 
this decay mode. If found, this would be another clear signal of this particular type of new 
physics. 

Z-mediated FCNC's can also contribute to the decays b — > at the one-loop level. 
However, the calculations are highly model-dependent, since it is necessary to include the 
new vector-singlet quark(s) with which the SM charge —1/3 quarks mix. The authors of 
Ref. [73] considered the case of a single vector-singlet quark, and included the contribution 
of the Higgs boson as well. For \U qb \ < 0.0017 they found that the contribution to b — > 
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was unimportant, but that the branching ratio for b — > d'j could be changed significantly, 
depending on the values of the masses of the exotic quark and the Higgs boson. Since 
there is quite a bit of model dependence, we do not consider this to be a clean signal of 
new physics. 

In summary, if Z-mediated FCNC's contribute significantly to B®-B® mixing, they will 
also lead to large effects in a variety of penguin decays. The present experimental upper 
limit on the FCNC parameters is \U q b\ < 0.0017 (q = d, s). This value for the parameters 
leads to unmistakable effects in b — > ql + l~ , B® — > Z + Z~, and electroweak penguins. In 
addition, Z-mediated FCNC's may lead to direct CP violation in decay modes such as 
B + -> n+n°. 



3.4) Multi-Higgs-doublet models 



(i) Natural Flavor Conservation 

In these models the new contributions to rare B decays come about through ampli- 
tudes in which charged-Higgs exchange replaces the SM W exchange. As noted in Eq. (28) 
and the surrounding discussion, the measurement of b — > s"f already excludes a region of 
parameter space with low M H + and high \Y\. 

In spite of this, charged-Higgs exchange may have significant effects on other rare B 
decays. The processes B® — > /+/", q = d, s were studied in multi-Higgs models in Ref. [74]. 
Neglecting small contributions from neutral pseudoscalar Higgs (which are proportional to 
m^/Mp, where M P is the pseudoscalar mass) [75], one finds the following expression for 
the ratio of charged-Higgs and SM amplitudes: 
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(63) 



where 
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(64) 



Since B(y t ) is negative, the two amplitudes add up constructively. 

To illustrate the effect, let us consider the case M H + — 400 GeV, Y = 3 discussed in 
Sec. 2.3(i). For these parameters, the ratio of Eq. (63) becomes 2.5. Thus, the B® — * l + l~ 
rates are expected to be an order of magnitude larger than in the SM, while charged-Higgs 
exchange dominates B^-B^ mixing. 
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The decays b — > ql + l~ in multi-Higgs-doublet models are more complicated than 
B® — > Z + Z~, since there are more operators which can contribute to this process. We 
refer to Ref. [76] for the details of the computation, but the conclusions are as follows. 
For those values of \Y\ and M H + for which B®-B® mixing is dominated by charged-Higgs 
exchange (i.e. large \Y\), the decays b — > ql + l~ can be enhanced by a factor of about 2. 
For b — > dl + l~ this is within the error of the SM prediction, but it is a significant effect 
for b — > s l + l~ . 

Finally, we turn to hadronic penguin decays. Model-dependent studies [77] show 
that, once the constraint from b — > is taken into account, the effect of charged-Higgs 
exchange can change the SM predictions by no more than a few tens of percent. Since 
these predictions suffer from large hadronic uncertainties, the rates of these processes 
cannot signal charged-Higgs effects. 

(ii) Flavor Changing Neutral Scalar s 

In models without NFC, one can have flavor-changing processes mediated by neu- 
tral scalars (FCNS). The flavor-changing couplings between quarks of flavor i and j are 
parametrized as (m,i/v)F i j('j 5 /2) (rrii > mj), while the flavor-conserving couplings are 
rrif/v. Thus there are also contributions to penguin decays due to neutral scalar exchange. 
Consider the decay b — > qff, where / can be a quark or a lepton. The rate for this decay, 
due to FCNS alone, is 

This rate is clearly maximized when the fermion / is as massive as possible. Consider then 
the decay b — > qr + r~. For M H o = 100 GeV, Fdb = 0.02 is essentially the largest value 
possible due to constraints from B^-B® mixing. On the other hand, F s i, has no similar 
constraint, so we take F s b = ^m s /mi, = 0.17. Then 

m^^r-)u K . ={^*;°:I ; p; (66) 

The branching ratio for b — > dr + r~ is only about 1/3 of the SM expectation, well within 
the errors of the prediction. This decay can therefore not be used to find effects of neutral 
scalars. On the other hand, the decay b — > st + t~ could be significantly affected by 
flavor-changing neutral scalars, since its new-physics branching ratio is of the same order 
as the SM prediction. However, note that we have selected almost maximal values for 
the new-physics parameters. If F s t, is smaller, or M H o larger, than the values we have 
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chosen, the new-physics contribution to b — > st + t~ would then diminish considerably, 
although there could still be important effects in B®-B® mixing. Note also that, due to the 
mass suppression, the FCNS contribution to decays involving lighter fermions is completely 
negligible. 

The one decay in which the mass suppression is not obviously a disadvantage is — > 
Z + Z~, q = d,s, since such a suppression is present even in the SM. We find 

T(B° q -> l + r)\ FCNS = -^-G 2 F M Bq f B mf \F qb \ 2 {^j . (67) 
Again taking M H o = 100 GeV, F sb = 0.17, and F db = 0.02, this gives 
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(68) 



For B® decays, the new-physics effects are within the errors of the SM prediction. For 
decays, the branching ratios due to flavor-changing neutral scalars are a factor of 2-3 times 
larger than in the SM. Since there are uncertainties in the SM prediction, and since we 
have taken optimal values for the new-physics parameters, this can only be considered a 
marginal signal of new physics. 

Therefore, in models with flavor-changing processes mediated by neutral scalars, there 
are no "smoking gun" signals in penguin decays. For maximal values of the new-physics 
parameters, there may be enhancements in the branching ratios of b — > st + t~ and S° — > 
/ + /~. However, for other values of these parameters there will be no significant effects in 
these and other penguin decays, even though there may still be important contributions 
to -Bg-Sq mixing. 

3.5) Supersymmetry 

The parameter space of supersymmetric models is quite complex, so that definite 
predictions of effects in penguin decays are hard to obtain. For example, consider the 
decay b — > S7 in the MSSM [35]. In addition to the charged-Higgs effects, which always 
increase the rate [see Sec. 2.3 (i)}, there are additional SUSY contributions from charginos 
+ up-type squarks or gluinos + down-type squarks in the loop. In certain regions of 
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parameter space the net effect is to cancel the contribution of the H ± , resulting in a 
branching ratio at or below the SM value, while in other regions the branching ratio is 
always greater than in the SM. Thus, the experimental measurement of b — > S7 does not 
constrain the SUSY parameter space in a simple way. 

The decay b — > sl + l~ has recently been analysed in supersymmetric models, taking 
into account the constraint from b — > S7 [78]. In the MSSM, it is found that the branching 
ratios for b — > se + e~ and b — > sfi + /j~ can be changed by at most 23% and 12%, respec- 
tively. These deviations are within the errors of the SM predictions, so these decay modes 
cannot be used to detect supersymmetry. However, the authors of Ref. [78] also study the 
C-odd lepton-antilepton energy asymmetry: 

A = N(Ei- >E l+ )-N(E l+ >E l -) 

N(E t - > E l+ ) + N(E l+ > Ei-) ' 1 > 

in which N(E[- > E t +) denotes the number of decays where the l~ is more energetic in 
the B meson rest frame than the Z + . They find that this asymmetry can be affected by up 
to 70% for b — > se + e~ and 48% for b — > s fi + \i~ . Furthermore, these sizeable deviations 
occur in a large region of SUSY parameter space. Thus, this asymmetry is an excellent 
place to look for effects of supersymmetry. Unfortunately, it is not clear how that region of 
parameter space which leads to large deviations in this asymmetry is correlated with that 
region of parameter space in which there are significant contributions to B®-B® mixing. 

Ref. [78] also examines b — > sl + l~ in a certain class of nonminimal SUSY models. In 
this case, the effects can be huge: the branching ratios for b — > se + e~ and b — > s fi + yT 
can be doubled, and the asymmetries enhanced by a factor of 3. 

For the decays B® — > and b — > qq'q' , similar studies have not yet been carried out 
in the context of supersymmetric models. However, since the branching ratio for b — > s 
is not substantially affected in the MSSM, this suggests that the SUSY contributions to 
B® — > and b — > qq'q', which are similar decays, will also not change the branching 
ratios significantly. On the other hand, in nonminimal SUSY models, the branching ratios 
for these decays may receive important corrections. 

4. Summary 

The phase information of the CKM matrix, which is the SM explanation of CP viola- 
tion, is represented by the unitarity triangle. At present, our knowledge of this triangle is 
rather poor - only the sides have been measured directly (IV^/V^I is probed in charmless 
B decays, and within the SM |Vtd/U c b| can be extracted from B^-B^ mixing), but these 
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measurements suffer from large theoretical uncertainties. In the near future, the angles of 
the unitarity triangle will be extracted from CP asymmetries in S-meson decays. Through 
the measurements of the CP angles a, f3 and 7, it will be possible to test the consistency 
of this description. There are three distinct ways in which the presence of new physics 
might be revealed: 

• (1): The relation a + P + 'j = n is violated. (Note that this relation can be tested only if 7 
is measured in CP asymmetries involving B® decays. If 7 is measured via B± _> D CP K± 
this relation will hold even in the presence of most types of new physics.) 

• (2): Although a + f3 + 7 = 7r, one finds values for the CP phases which are outside of 
the SM predictions. 

• (3): The CP angles measured are consistent with the SM predictions, and add up to 
180°, but are inconsistent with the measurements of the sides of the unitarity triangle. 

If any of these discrepancies is found, we will want to know what kind of new physics 
is involved. The principal way in which new physics can affect the CP asymmetries is 
through new contributions to B®-B® mixing. By performing a model-by-model analysis 
of physics beyond the SM, it is possible to ascertain which types of new physics might be 
responsible for the discrepancy. This analysis allows us to separate new-physics models 
into two types: (i) those in which the phase of B®-B® mixing is changed, in which case 
each of items (l)-(3) may occur, and (ii) those in which the phase is unchanged, in which 
case only (3) is possible. However, this analysis does not tell us how to distinguish among 
models of a given type. It is this question which we have attempted to address in this 
paper. 

Our main observation is quite simple. Any new physics which contributes to B^-B^ 
or B®-B® mixing will necessarily contribute to the rare flavor-changing penguin decays 
b — > dX and b — > sX. In some cases, the values of the new-physics parameters which 
yield significant effects in B®-B® mixing will also lead to large deviations from the SM 
predictions for certain penguin decays. It is therefore possible to partially distinguish 
among different models of new physics by examining their predictions for these penguin 
decays. 

• (i) Four generations: The phase of B®-B® mixing can be changed due to the new box- 
diagram contributions with internal t' quarks. Although the new-physics parameter space 
is too large to make absolute predictions for branching ratios of penguin decays, we have 
shown that there are regions of parameter space in which both B CP asymmetries and B 
penguin decays are significantly affected. The particular example we chose found important 
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contributions to B^-B^ mixing, and roughly an order-of-magnitude enhancement of the 
branching ratios for both exclusive b — > d electroweak penguins and B® — > 

• (ii) Z -mediated flavor-changing neutral currents: The phase of B®-B® mixing can be 
altered due to the tree-level exchange of a Z with flavor-changing couplings. In fact, B®- 
B^ mixing can be dominated by this new physics. If these new contributions are important, 
there will also be unmistakable effects in b — > ql + l~ , B® — > / + /~, and electroweak penguins. 
In particular, the rates for b — > s (b — > d) penguin processes can be enhanced by as much 
as one (two) orders of magnitude. In addition, Z- mediated FCNC's may lead to direct CP 
violation in decay modes such as B + — > 7r + 7r°. On the other hand, if the branching ratios 
for b — > ql + l~ are found to be consistent with the SM, this will indicate that Z-mediated 
FCNC effects in B®-B® mixing are negligible. 

• (Hi) Multi-Higgs- doublet models with natural flavor conservation: There are new box- 
diagram contributions to B®-B® mixing with internal charged Higgs bosons, but the phase 
of this mixing is unchanged. (If one also has spontaneous CP violation in such models, 
this phase is zero, since the CKM matrix is real, and the unitarity triangle becomes a 
straight line.) For that region of parameter space in which B®-B® mixing is dominated 
by the charged-Higgs contribution, the branching ratios for B® — > and B — > X q l + l~ 
are enhanced by up to an order of magnitude or a factor of 2, respectively. However, 
when the box diagrams with internal and bosons are about equal in magnitude, 
charged-Higgs effects in penguin decays may not be sufficiently large to be detected, due 
to theoretical uncertainties in the SM rate calculations. 

• (iv) Multi-Higgs- doublet models with flavor-changing neutral scalars: The phase of B®- 
B® mixing can be changed due to the tree-level exchange of a neutral scalar with flavor- 
changing couplings. However, there are no significant effects in B penguin decays. This 
is due to the fact that the flavor-conserving coupling of the neutral scalar to a fermion is 
proportional to the fermion mass, and penguin decays all involve light fermions. 

• (v) Left-right symmetric models: There are no significant effects in the B system in these 
models. The only exception is the case where the right-handed CKM matrix is considerably 
fine-tuned, but we do not consider this possibility. 

• (vi) Minimal super symmetric models: There are many new contributions to B®-B® mix- 
ing, but all have the same phase as in the SM. A search of the parameter space reveals that 
SUSY contributions to b — > sl + l~ do not change its branching ratio significantly. This 
suggests that the branching ratios for the decays B® — > / + /~ and b — > qq'q' will also be 
relatively unaffected. However, SUSY can be detected by examining the lepton-antilepton 
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energy asymmetry. This asymmetry can be affected by up to 70% for b — > s e + e and 48% 
for b — > s fi + fx~. 

• (vii) Nonminimal super symmetric models: In nonminimal SUSY models with quark- 
squark alignment, the SUSY contributions to B®-B® mixing (and hence to B penguin 
decays) are generally very small (though models do exist in which Mf 2 c/SY '(£§) /Mf 2 M \B° d ) « 
0.15). In alternative nonminimal SUSY "models" it is simply assumed that all SUSY 
parameters take the maximum allowed values allowed by experiment. These models are 
not terribly predictive, due to the very large numbers of parameters. It is possible to find 
new contributions to the mixing with different phases than in the SM, and to arrange 
the many parameters such that the branching ratios of penguin decays are enhanced or 
suppressed. However, it is virtually impossible to analyse the effects of nonminimal SUSY 
on B®-B® mixing and B penguins in any systematic way. 

If some indication of new physics is found in the measurements of CP asymmetries, 
the above analysis may be used to distinguish different candidate models of physics beyond 
the SM. (In fact, in some cases it is likely that the new physics will be found first through 
measurements of rare B decays.) For example, suppose that new physics is found through a 
discrepancy of type (1) or (2). This would indicate that the new physics is probably either 
a fourth generation, Z-mediated FCNC's, or flavor-changing neutral scalars. Since each 
of these three models affects B penguin decays differently, they can be at least partially 
differentiated by a study of such decays. And if the new physics is found through a 
discrepancy of type (3), the new physics is likely to be either a multi-Higgs-doublet model 
with NFC, or the MSSM. In this case, it may be difficult to distinguish the two models 
of new physics since their effects on B penguin decays are similar. This is not surprising, 
since the MSSM contains two Higgs doublets. Still, there are signals, such as the lepton- 
antilepton energy asymmetry in b — > which can differentiate these two models. 

To sum up, most physics beyond the SM which can affect CP asymmetries in B decays 
will also contribute to rare, flavor-changing B decays. We have examined the effects of a 
number of models of new physics on both CP asymmetries and penguin decays. Although 
not all models of new physics have "smoking gun" signatures in these decays, we have shown 
that the measurements of CP asymmetries and rare penguin decays give complementary 
information, and both will be necessary if we hope to identify the new physics. 
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